Draft version January 16, 2012 

Preprint typeset using L^I^X style cmulatcapj v. 5/2/11 



(N 

o 

a 

(N 



< 

o 

43 

Oh- 
i < 

o 

CO 

03 



> 
OO 
(N 

o 

(N 



13 



NOMADS OF THE GALAXY 
Louis E. Strigari 1 , Matteo Barnabe 1 , Philip J. Marshall 2 , Roger D. Blandford 1 

Draft version January 16, 2012 

ABSTRACT 

We estimate that there may be up to ~ 10 5 compact objects in the mass range 10~ 8 — 10~ 2 Mq 
per main sequence star that are unbound to a host star in the Galaxy. We refer to these objects as 
nomads; in the literature a subset of these are sometimes called free-floating or rogue planets. Our 
estimate for the number of Galactic nomads is consistent with a smooth extrapolation of the mass 
function of unbound objects above the Jupiter-mass scale, the stellar mass density limit, and the 
metallicity of the interstellar medium. We analyze the prospects for detecting nomads via Galactic 
microlensing. The Wide-Field Infrared Survey Telescope ( WFIRST) will measure the number of 
nomads per main sequence star greater than the mass of Jupiter to ~ 13%, and the corresponding 
number greater than the mass of Mars to ~ 25%. All-sky surveys such as Gaia and LSST can identify 
nomads greater than about the mass of Jupiter. We suggest a dedicated drift scanning telescope that 
covers approximately 100 square degrees in the Southern hemisphere could identify nomads as small 
as 1O~ 8 M via microlensing of bright stars with characteristic lightcurve timescales of a few seconds. 
Subject headings: gravitational lensing - planets and satellites: general - Galaxy: general 



1. INTRODUCTION 

The recent years have witnessed a rapid rise in the 
number of known planetary mass objects, < O.OlAf , 
in the Galaxy. Searches for exoplanets from radial ve- 
locities find that ~ 30 — 50% of GK dwarf stars have 
planets greater than the mass of Neptune within peri- 
ods < 50 days (| Wolfgang fe Laughlhfeoill) . Transiting 



searches find that ~ 15% of main sequence dwarfs are or- 
bited by short-period planets at less than four Earth-radii 
(jBorucki et al.ll2010f ). Direct imaging and microlensing 
have a lso now start ed to uncover planets bound to host 
stars (IGaudill2010h . 

Much less is known about the population of < 0.01M Q 
objects that are not bound to a host star. Several can- 
didate unbound objects have been imaged in star clus- 
ters with mass possibl y as small as a few times that o f 
Jupiter, 10~ 3 M(T) (e.g . I Javawardhana fe Ivanovl (|2006h : 
iCaballero etTal] (|2007t) : iBihain et all (|2009t) L However, 
the origin of these objects is uncertain; they may 
have f ormed dire ctly in the collapse of the molecular 
cloud (Recs 1976), or have been ejected from their birth- 
lace aroun d a host star via a dynamical interaction 
Bossll2000f) . Free-floating objects at the Jupiter mass 
and below have been difficult to find via these methods. 

Microlensing, however, does provide a way — and 
perhaps the only way — to detect objects below the 
deuterium-burning mass limit that are not bound to a 
host star. In a recent sur vey of the Galac tic bulge, 
the MOA-II collaboration (jSumi et al.1 120 111 ) reported 
the discovery of planetary-mass objects either very dis- 
tant from their host star (~ 100 AU) or unbound from 
a host star entirely. This detection was obtained from 
analysis of the timescale distribution of the microlensing 
events, which showed a statistically-significant excess of 
events with timescale < 2 days as compared to a stan- 

1 Kavli Institute for Particle Astrophysics and Cosmology, 
Stanford University, Stanford, CA 94305 USA 

2 Department of Physics, University of Oxford, Keble Road, 
Oxford, OX1 3RH, UK 



dard Galactic model with a stellar mass function cut-off 
at the low mass end of the brown dwarf regime. The 
mass function of this new population of objects can be 
described (for illustration) by a <5-function with a best fit 
near the Jupiter mass. These results tell the surprising 
story that objects greater than about the mass of Jupiter 
are approximately twice as numerous as both main se- 
quence stars and planets bound to host stars. 

Though their existence is established, the origin of 
these unbound objects is far from clear. Do they form a 
continuation of the low end brown dwarf mass function 
near the deuterium burning mass limit, or did they form 
as a distinct population of objects ejected from their orig- 
inal host stars? Because of their uncertain origin and 
their present status, we prefer to refer to objects with 
mass < O.OIMq that are not bound to a host star as 
Nomads; in the literature they have been also referred 
to as rogue or free-floating planets. The name "nomad" 
is invoked to include that allusion that there may be 
an accompanyi ng "flock," either in the form of a sys- 
tem of moons (jDebes fe Sigurdssod[2007l ) or in its own 
ecosystem. Though an interstellar object might seem an 
especially inhospitable habitat, if one allows for inter- 
nal radioactive or tectonic heating and the development 
of a thick atmosphere effective at trappi ng infrared heat 
(|Stevensonlli999l : lAbbot fc Switzerll2011[ ) , and recognizes 
that most life on Earth is bacterial and highly adaptive, 
then the idea that interstellar (and, given the prevalence 
of debris from major galaxy mergers, intergalactic) space 
is a vast ecosystem, exchanging mass through chips from 
rare direct collisions, is intriguing with obvious implica- 
tions for the instigation of life on earth. 

Understanding the bounds on the nomad population, 
and the prospects for detecting them with microlensing 
surveys, is the focus of this paper. In particular, what is 
the number and mass density of nomads in the Galaxy? 
What are the bounds on the minimum mass of a nomad, 
and what is the lightest detectable nomad? How well can 
we measure the nomad mass function, and how does this 
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compare to the low mass brown dwarf mass function? 
And can we independently measure the mass function of 
nomads in the bulge and in the disk? 

We will show that a dedicated space-based survey of 
the inner Galaxy, such as the proposed Wide-Field In- 
frared Survey Telescope (WFIRSlE), will measure the 
number of nomads per main sequence star greater than 
the mass of Jupiter to ~ 13%, and the corresponding 
number greater than the mass of Mars to ~ 25%. We 
also show that large scale Galaxy surveys, in particu- 
lar the GaiaQ mission and the Large Synoptic Survey 
Telescope (L5570), will be sensitive to nomads greater 
than about the mass of Jupiter without changing their 
proposed observing plan. Further, the back-to-back 15 
second exposures in the planned design of LSST will al- 
low for limits at least to be placed on nomads near the 
mass of Pluto, < 10 _8 Mq. As an extension, we suggest 
that a dedicated drift scanning telescope could identify 
nomads as small as 1O~ 8 M0 via microlensing of bright 
stars with characteristic lightcurve timescales of a few 
seconds. 

This paper is organized as follows. In Section[2]we esti- 
mate the number of nomads in the Galaxy. In Section [3] 
we calculate the event rate of nomads in microlensing 
surveys. In Section 3] and Section [5] we outline methods 
for measuring the nomad population and simulating de- 
tection efficiencies. In Section[6]we present the results of 
these projections. In Section [7] we postulate a survey for 
short timescalc nomads via a drift scanning telescope. In 
Section we summarize our conclusions. 

2. THE NOMADIC POPULATION 

We begin by setting up the model for the nomad pop- 
ulation. Objects with mass < 10~ 2 Mq are believed 
to originate from two distinct processes. Between the 
Jupiter mass and the deuterium-burning mass, many 
of these objects may form similar to stars by gravita- 
tional fragmentation. Below Jupiter masses, they likely 
are born in protoplanetary disks and dynamically-ejected 
during the evolution of the system. It is unknown from 
a theoretical perspective whether there is a smooth con- 
tinuation of the mass function at the dividing mass that 
separates these populations. 

In light of these uncertainties, we choose a simple 
broken power-law model for the nomad mass function, 
dN/dM oc M~ a , which is a smooth continuation of the 
mass function at higher masses, 

a nm for 1(T 8 M/Mq < 0.01 

= J a b d for 0.01 < M/Mq < 0.08 

a 2 for 0.08 < M/M Q < 0.70 

ai for 0.70 < M/M Q . 

In addition to these power law slopes, we define the min- 
imum cut-off in the nomad mass function as m m i n , the 
number of objects in the nomad mass regime as iV nm , and 
the number of main sequence stars in the mass regime of 
0.08 — IMq as -/Vms- From the latter two quantities we 
define the ratio f3 = N nm /Nus- 

For the above parameterization, empirical bounds over 
the entire nomad range may be motivated from several 
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considerations. First, the mass function of the lowest 
mass nomads that we consider may be estimated from 
bounds on the population of Kuipcr Belt objects (KBOs). 
We start fro m the distribution of d iameters of KBOs de- 
termined in iBernstein et al.1 <|2004f ) . At the high diam- 
eter end, D > 100 km, the KBO distribution scales as 
dN/dD cx D- 4 . Below the break radius of D ~ 100 km, 
where collisional effects are believed to be important, the 
KBO distribution flattens, dN/dD oc const. Assuming 
that the bodies of the outer Solar System have approx- 
imately constant mass density of ~ 1 g cm -3 , the mass 
function scales as dN/dM oc m~ 2 above 10 _12 Mq, while 
below the mass function scales as dN/dM oc vnT 2 ! 3 . 

At the highest mass end, corresponding to approxi- 
mately several times the mass of Jupiter, there is evi- 
dence that nomadic objects in open clusters constitute a 
smooth continuation of the bro wn dwarf mass functio n 
at higher masses with a = 0.6 (|Caballero et al.ll2007|) . 
Further there may be a turn over in the mass fun ction be- 
low ~ 6 times Jupiter mass (jBihain et all 120091 ) . though 
these results are subject to systematic uncertainties on 
the masses of the objects and the small number of objects 
known. 

In comparison to these results from direct imaging, 
microlensing observations more strongly constrain the 
nomadic mass function, in particula r at the hig h mass 
end. The microlensing results from iSumi 
find that the equivalent best-fitting slopes and one-sigma 
uncertainties are a nm = 1.3^04 an d Q!bd — 0.481°; 2 ^. 
In lSumi et al.1 ()2011h the minimum mass was taken to be 
w m i n = 1O _5 M0, though given the cadence of the sur- 
vey they are insensitive to values of m m i n at this mass 
scale and below. Taking ro m i n = 1O _5 M0 and these 
best-fitting slopes implies /3 ~ 5. Extrapolation down 
to below Earth mass scales, m m - m = 1O _6 M0, yields 
(3 ~ 10, and further extrapolation down to m mm = 
1O _8 M0 yields (3 ~ 60. Intriguingly for a continuous 
power law extrapolation down to ~ 1O~ 15 M0, the num- 
ber of nomads per star approaches the bound on the 
abun dance of interstellar comets (jFrancis! 120051 : Uural 
120 lit ), and the corresponding nomad mass density is 
~ 1% o f the oxygen mass density in th e interstellar 
medium (|Baumgartner fc~M ushotzkv 2006T). 

Assuming the above parameterization of the mass func- 
tion, a nm is negatively correlated with from mi- 
crolensing observations. For example, for abd = E the 
95% c.l. lower limit on the nomad slope is a nm = 0.5. 
For this combination of slopes, extrapolating down to 
n^min = 10~ 8 Ajf Q impl ies that (3 J> 1 . On the other 
hand for a hd = 0, the ISumi etaLI (|2011l ) 95% c.l. up- 
per limit on the slope is a nm = 2. In this case assum- 
ing m m i n = 1O _5 M0 implies /3 ~ 50. Extrapolation 
down to m m in = 10 _6 Mq implies an order of magni- 
tude increase in ft ~ 700, while extrapolation down to 
77i min = 10 _8 M o implies f3 ~ 10 5 . 

The above estimates indicate that, when fixing to 
the measured abundance of nomads at > 1O~ 3 M and 
smoothly extrapolating to lower masses, there is several 
orders of magnitude uncertainty on the nomad abun- 
dance. For an appropriately large ratio of the number 
of nomads to main sequence stars, the nomad mass func- 
tion is constrained by the limits on the number of com- 
pact objects in the Galactic disk and halo. For masses 
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Figure 1. Number of nomads greater than a given mass scale, 
N n m(> M), relative to the number of main sequence stars, A^ms- 
Three different slopes for the nomad mass function are labeled, 
c»nm = 2, 1.3,0.5. The upper black curve is determined assuming 
that objects at that mass scale have a density of po = O.IMq pc~~ ' ! . 

> 10 _7 M Q , the abundance of compact objects in the 
halo is constrained by uppe r limits on MACHO dark mat- 
ter (jTisserand et alJ l2007f ). These bounds indicate that 
compact objects of mass > 1O _7 M comprise < 10% 
of the Galactic dark matter halo. For more local mea- 
surements, a less stringent limit arises from the stellar 
mass density of th e Galactic disk, which we take to be 
po = 0.1M@ pc~ 3 (|Holmberg fc Flvndl2000l) . As an ex- 
ample, if we assume that the entire local mass distribu- 
tion of the disk is comprised of objects at the mass scale 
1O _8 M0, the bound po = 0.1M Q pc -3 corresponds to an 
upper limit of 10 6 compact objects per main sequence 
star. Masses of compact objects at these scales and be- 
low may be probed by future short cadence microlensing 
obser vations with the Kepler satellite (see iGriest et al.l 
120111 and the discussion below). 

Predictions for the number of nomads, in comparison 
to constraints on the local mass density, are summarized 
in Figure [TJ This shows the number of nomads greater 
than a given mass scale, N(> M), relative to the num- 
ber of main sequence stars, Nms- Four different slopes for 
the nomad mass function are labeled, a nm = 2,1.3,0.5, 
which have corresponding values for the slope of the 
brown dwarf mass function of abd = 0,0.5, 1. The up- 
per limit, indicated as the diagonal line, is determined 
assuming that objects at that mass scale have a density 
of po = O.IMq pc -3 . For a nm > 2, the total mass of 
the nomad population is dominated by the lowest mass 
objects. 

3. EVENT RATES 

In this section we calculate the microlensing event rates 
from the nomad population. We begin by establishing 
the definitions and the Galactic model parameters, and 
then use this model to predict the timescale distribution 
of events and the integrated number of events detectable. 

3.1. Definitions 



We employ standard microlensing formalism. The dis- 
tance to the source star is D$, the distance to the lens is 
Dl, and the mass of the lens is M. The Einstein radius 
is i?| = (4G/c 2 )MD L (D s - D L )/D S , and the Einstein 
crossing timescale is tE = Re I v. The amplification of 
a source star is A{u) = (u 2 + 2)/(u\/u 2 + 4), where u is 
the projected separation of the lens and source in units 
of the Einstein radius. 

To calcu late microlen s ing ev ent rates we use Galactic 
model 1 of iRahal et al.l (|2009() . This is characterized by 
an exponential thin disk with a scale length of Rd = 3.5 
kpc and scaleheight of Zh = 0.325 kpc, 



p d {R) = p exp[-(i?- R )/Rd - \z\/zh] 



(1) 



where R = 8.5 kpc, %% = 0.35 kpc. We us e the bulge 
density distribution from IDwek et al.l (|1995l ). The lens- 
source transverse velocity distribution, f(vi,Vb) is mod- 
eled as in IHan fc Gouldl (fl996l) . with v = y/vf+v 2 , 
where we indicate as vi and i>&, respectively, the velocity 
along the galactic longitude and latitude coordinates. 

We determine the total event rate distribution by 
breaking the lenses-sources into the disk-bulge, bulge- 
bulge, and disk-disk components. We consider two differ- 
ent targets of source stars. First, bulge stars in the direc- 
tion of Baade's window, (b, I) = (—3.9°, 1°), and second, 
sources distributed over all-sky. For the bulge observa- 
tions we can compar e to observationa l determinations 
of th e optical depth dSumi et al.l 120031; iPopowski et al.l 
120051: iSumi et al.l 120061: lHamadache et all 120061) and "to 



the theoretical optical de pth calculations (jHan fc Gouldl 
120031 : iWood fc Maoll2005] ) by smoothly extrapolating the 
rates from the nomad mass regime to the mass regime of 
main sequence stars and remnants. 

3.2. Finite source effects 

Since we extrapolate the nomad mass function down 
to low mass scale, it is important to properly account 
for finite source effects in the microlensing events. More 
specifically, we need to estimate by how much the peak 
amplification of an event is reduced when Re is of order 
the projected radius of the source star. To estimate finite 
source effects we take the source stars to have uniform 
surface brightness, and for a given projected lens-source 
separation we estimate the amplification as 



A is (u) 



2tt 



dej) pdpA(\J u 2 + p 2 — 2upcos4>), 



(2) 

where p* = Rq/Re- For typical lens and source dis- 
tances, Dl ~ 5 kpc and D$ — 8 kpc, and assuming 
a lens mass M = 10 _8 M Q , the peak amplification is 
Af s ~ 1.1. Though this is less than the standard point 
lens-point mass amplification by ~ 15%, surveys that we 
consider below will still be sensitive to brightness fluctu- 
ations of this magnitude. Extrapolating further down to 
M = 1O _9 M0, the peak amplification is only A{ s ~ 1.01. 
Though extraction of events at this brightness may still 
be detectable, to provide conservative estimates we re- 
strict our analysis to lens masses > 1O _8 M0. 

3.3. Bulge event rate 

The microlensing event rate per source star in a di- 
rection (6, 1) is given by an integral over the lens- 
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source transverse velocity distribution, the lens density 
distributions p E and the mass function (jGriestl 119911: 
iKiraga fc PaczvnskilH99l . 

—jj- 2 — =«t J dD L J dvidv b vf(v b ,vi)8(t E - R E /v) 



(3) 



x / dM((M)R EPL (l,b,D L ). 



The mass function £(M) is normalized to the mean mass 
of the lens population. The optical depth is r ~ 7tT/2, 
where T is the integral of the event rate distribution over 
all t E . In Eq. 1 u T = (A 2 /VA 2 - 1 - l) 1 / 2 , with A = 
1.34 corresponding to ut = 1. This corresponds to the 
event rate for source stars within a circular area of one 
Einstein radius of the lens star. This is a conservative 
criteria that is appropriate for our analysis; event rates 
are increased for ut > 1 when allowing for A < 1.34. 

Figure [2] shows the event rate distributions in the di- 
rection of the bulge, with each panel corresponding to 
a different assumption for the slope of the nomad mass 
function. In all panels we take ct\ = 2.0 and U2 = 1.3 for 
the main sequence stellar mass function. Within each 
of the panels there are three different assumptions for 
Timin- The three curves in all of the panels represent the 
sum of the event rate distribution from bulge and disk 
lenses. For all curves in the middle and right panels, the 
mean timescale of a microlensing event is (t E ) ~ 50 days. 
However, for the curves in the left panel, (t E ) depends 
strongly on m m i n because of the steep power law to low 
masses. Specifically for m min = (10~ 2 , 10~ 5 , 1O _8 )M , 
the respective mean timescales are (t E ) ~ (50,35,3) 
days. For all curves, the optical depths are ~ 1.5 x 10~ 6 , 
consistent with the theoretical calculations above and 
the observations. The inclusion of the nomad popula- 
tion does not affect the total optical depth because this 
quantity is independent of the mean mass of the lens 
population. 

3.4. All-sky event rate 

We now move on to examining the all-sky event rate 
distribution. In addition to the ingredients input into 
Eq. [3j here we require two additional pieces of informa- 
tion: the luminosity function of sources, <f>(m), where 
m is the source apparent magnitude, and the radial dis- 
tribution of sources ng(r). For the former we use the 
solar nei ghborhood y-band luminosi ty function as com- 
piled in iBinnev fc Merrifieldl (|1998f ). and the I^-band 
du st extinction model for th e Galaxy as parameterized 
in iBelokurov fc Evansl (|2002h . For the latter, we scale 
the disk density profile by the local density p<j, i.e. 
n s = Pd/po- Note that here we exclude bulge sources 
because they only have a few percent contribution to the 
all-sky microlensing event rate. 

With the definitions above, the total microlensing 
event rate brighter than a limiting magnitude, mi; m , is 

/" m iim r°° 
r(< TOiim)=MT / (f)(m)dm dD s D 2 s n s (l,b,D L ) 
Jo Jo 
rDs r 
x J dD L J dvidv b vf(v b ,vi) 



x / dMC(M) f R E p L (l,b,D L ) 

J m m ; n J 



(4) 



Figure [3] shows the integrated all-sky event rate as a 
function of the limiting magnitude, for the same sets of 
nomad and brown dwarf mass function parameters that 
are shown in Figure [2] Here we have included only the 
event rate for 30 minutes < t E < 1 day. The lower bound 
for this timescale distribution is motivated by considering 
the mean timescale for an object of mass 10 _8 Mq, while 
the upper bound is motivated from Fig. [2] which shows 
that events from objects with mass < 1O -2 M0 predomi- 
nantly have timescales < 1 day. We will further motivate 
the lower cut-off of 30 minutes when we discuss analy- 
sis of observational prospects below. For each curve, the 
value of TO m i n is indicated. As Figure |3] shows, there is 
~ 4 orders of magnitude uncertainty in the predicted no- 
mad event rate brighter than 20th magnitude. For the 
most shallow allowable nomad mass function, a nm = 0.5, 
the event rate in this range of timescales is ~ 0.2 per year, 
while for the steepest allowable mass function, a nm = 2, 
the event rate is > 10 3 per year. 

We note that, when restricting to lens masses > 
10 _2 Mq, the event rates determined in Fig. [3] are con- 
sistent with prior estimates of ~ few per year for V < 15 
(|Nemiroffl II 9981 : iHanl 120081 ) . Including the entire pop- 
ulation of nomads, stars, and remnants, in fact we esti- 
mate ~ 2500 photometric microlensing events for sources 
greater than 20th magnitude. Again the vast majority of 
these events are from disk sources from the high density 
region towards the Galactic center, with a few percent 
contribution from bulge sources. The challenge for fu- 
ture observations will clearly be to achieve the appropri- 
ate efficiency to extract these short timescale events. 

4. FORECAST METHODOLOGY 

The results from the section above provide an estimate 
of the nomad event rate, independent of the survey spec- 
ification. In this section, we use the above predictions 
to estimate how well the nomad population can be mea- 
sured, given some basic input variables for a survey. 

As a general strategy, we would like to determine the 
constraints on a nm , abd and m m i n likely to be available 
from surveys of varying cadence, exposure, and sky cov- 
erage. Here we define the exposure in a standard manner 
as the number of stars monitored, JV*, during an observa- 
tional time period, T Q b s - For the given exposure and the 
Galactic model discussed above, we take the data as the 
observed timescale distributions for a set of microlensing 
events. We assume that there are n bins distributed over 
the range of observed t E . The minimum and maximum 
detectable timescale for a survey is set by the detection 
efficiency, which we estimate below for surveys of differ- 
ent cadence and exposure. 

We denote T t as the mean event rate in the i th t E 
bin for a specified exposure, where I\ is a function of 
the model parameters a nm , «bd and m mm . Our goal 
is to estimate how well we can measure these parame- 
ters for an observed event rate distribution. All of the 
other parameters, such as the local stellar density, the 
disk scale length and scale height, the bulge mass distri- 
bution, and stellar mass function in the regime of main 
sequence stars and above are fixed to their fiducial val- 
ues. We do this primarily for simplicity in order to effec- 
tively isolate the impact of the nomad population. We 
assume that the probability for obtaining N t events in 
the i th timescale bin follows a Poisson distribution with 
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Figure 2. The event rate timescale distribution for several mass functions and cut-off masses. In all panels the solid, dotted, and dashed 
curves assume rn m \ n = 10~ 2 , 10 -5 , 1O~ 8 M0, respectively. The slopes of the mass function in the planetary mass regime (a n m) are indicated 



curves assume m mln 
in each of the panels. 



a mean \x % = T ^ s NnT l . For the assumptions above, we 
can define the elements of the inverse covariance matrix 



as 



T ohs N* dT l dT l 

r, 86a de b ' 



(5) 



where the indices a and b represent the model parame- 
ters, which in our case are a nm , abd, and m m i n . From 
Eq. [5j the one-sigma uncertainty on parameter a is F"^ , 
evaluated at the fiducial values for the parameters. To 
evaluate Eq. [5] we choose the number of bins n to be 
equally spaced in log intervals. The main constraint on 
the bin size will be to ensure that they are wide enough to 
accommodate a 50% uncertainty in the reconstructed i^. 

5. DETECTION EFFICIENCY 

In the analysis above, we assumed 100% efficiency 
when detecting nomads over the entire range of event 
timescales. In order to obtain a more realistic event rate 
for a specific survey, we must gain an understanding of 
how the detection efficiency scales as a function of event 
timescale. In this section, we describe the basic set-up 
for our efficiency simulations, and how they are adapted 
to specific surveys in the sections that follow. 

We estimate the detection efficiency via a basic pro- 
cedure for generating microlensing events. We begin by 
drawing source and lens objects from the appropriate 
disk or bulge density distribution. The relative trans- 
verse vejflcjtv_js_then_drawn from the velocity distribu- 
tion (jHan &: Gould fl996l ). We draw the impact parame- 
ter for the source and lens randomly on a uniform inter- 
val out to the Einstein radius, and the peak timescale of 
the event, to, uniformly during the duration of a given 
survey, T obs . 

The above set of parameters, {Ds, Dl, v, to), along 
with the event timescale ts fully describe the microlens- 
ing event. For this set of parameters, we compute the 
amplification of the source as a function of time, which by 
definition peaks at to. The amplification is calculated at 
time-steps specified by the cadence of the survey. Moti- 
vated by the two different survey set-ups that we discuss 
below, we consider two different models for the survey 



cadence. First, we consider a uniform cadence model 
in which the number of time-steps is simply T b s /t ca .d, 
where i ca d is the cadence of the survey. Second, we con- 
sider a quasi-uniform cadence, in which there are a total 
of n e epochs for the survey, and n m measurements uni- 
formly spaced per epoch. As discussed below this is most 
relevant when discussing results for the Gaia survey. 

For each point on the lightcurve, the error is estimated 
from the expected photometric precision. Since the focus 
of our analysis is on bright microlensing events, we take 
the error to be uniform for all source stars, and character- 
istic of the survey that is considered. We then simulate a 
lightcurve point by sampling from a normal distribution 
centered on the true point with a variance given by the 
photometric error. The specific error assumed for each 
survey will be provided below. 

With the above procedure in place, it remains to quan- 
tify a criteria for detection for a microlensing event; we 
choose a relatively simple one that is appropriate for the 
scope of this work. For our primary analysis we de- 
mand that three consecutive points on the lightcurve are 
> 3<7 deviations from the mean baseline m agnitude of a 
star This has been used in previous studies (jGriest et al.l 
120 111) , and p rovides a good app roximation to the criteria 
discussed in lSumi et al.l (|2011[ ). The detection efficiency 
for an input timescale is then the ratio of the number of 
simulated events that pass the selection criteria to the 
total number of simulated events at the input timescale. 

6. PROJECTIONS AND CONSTRAINTS FOR SPECIFIC 
SURVEYS 

With the above ingredients in place, we now move on 
to discussing event rates and constraints for specific sur- 
veys. We begin by examining next-generation bulge sur- 
veys with WFIRST, and then move on to discuss forth- 
coming all sky surveys Gaia and LSST. We conclude by 
examining the detection prospects in the short term for 
the Kepler satellite. 

6.1. WFIRST 

We first consider the case of a dedicated survey to mon- 
itor the inner Galaxy region. This is similar in spirit to 
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Figure 3. Event rate for sources brighter than a limiting magnitude, for different assumed mass functions and cadences. Left panel assumes 
the nomad and brown dwarf mass function is described by a n m = 2 and ct b d = 0; middle panel assumes cr n m = 1.3 and c*bd = 0.48; right 
panel assumes o nm = 0.5 and ct^d = 1.0. In each panel, the event rate is integrated between timescales of 30 minutes < tE < 1 day, and 
the value of m m ; n for each curve is indicated in the upper right panel. 



the modern MOA, OGLE, and EROS surveys, and to a 
larger scale, space -based extension s uch as the proposed 
WFIRST mission (|Green et alj|2011l ). For the former set 
of surveys, we can directly use their published detection 
efficiencies to predict the event rates and model the er- 
ror distributions, while for a WFIRST type mission this 
requires simulating events as described above. 

For WFIRST, we use a cadence of 15 minutes, a to- 
tal exposure time of 1 year, and photometric errors of 
0.1%, which will be achievable down to J = 20.5. Using 
the above model, at tE ~ 0.03 days we find a detec- 
tion efficiency of ~ 50%. This high efficiency at short 
timcscale is primarily driven by the order of magnitude 
increase in the photometric precision relative to mod- 
ern microlensing experiments. We note that if we as- 
sume the MOA-II cadence and photometric uncertainty 
in their high cadence fields, which we approximate as 
~ 30 observations per night at ~ 15 minute cadence, at 
t B = (0.5, 1, 10) days, we find efficiencies of (10, 20, 40)%, 
which prov ides a good approx imation to the efficiencies 
reported in iSumi et~aLl (|20P1 . 

In Figure 0] we show the resulting one- and two-sigma 
uncertainties on the combination a nm -ahd, for modern 
and for future dedicated surveys. Here we have assumed 
fiducial values of a^d = 0.48 and a nm = 1.3, though 
we generally find that our results are insensitive to the 
specific value for these quantities. In plotting the unfilled 
contours in the left panel, we have assumed an exposure 
and detection efficiency similar to the MOA-II analysis, 
which provides a total of ~ 500 events for 2 years of 
observations of 50 million stars; we have assumed n = 20 
bins distributed uniformly in log between timescales of 
1-200 days. In this case the errors from our model are 
in good agreements with the one-sigma unc ertainties on 
ctbd and a nm presented in ISumi et all (|2011f) , with slight 
departures due to the non-gaussian behavior in the tails 
of the results from the later. In this case the one-sigma 
errors arc c Qbd ~ 0.30 and <r Q , nm ~ 0.40. 

The filled set of contours in the left panel of Figure |4] 
show the projected constraints assuming a cadence of 
15 minutes and 2 x 10 8 monitored stars for one year. 



This cadence and exposure is mot ivated by the prelim - 
inary specifications for WFIRST ([Bennett et all 12010( 1 . 
To provide the most optimistic predictions, and as moti- 
vated by the photometric precision and the simulations 
described above, here we have assumed a 100% detec- 
tion efficiency at all tE > 0.04 days. In this case, the 
one-sigma uncertainties are reduced to cr antn = 0.03 and 
fr Qbd ~ 0.05, representing nearly an order of magnitude 
improvement relative to the modern constraints. If we 
assume m m ; n = 1O _3 M , this corresponds to a mea- 
surement of to ~ 13% precision, and for m m ; n = 
3 x 10 _7 Af Q we have a measurement of (3 to ~ 25% pre- 
cision. 

For comparison to the bulge results, in the right panel 
of Figure Q] we show the resulting constraints for disk ob- 
servations towards (—2.4°, 331°). This direct ion is specif- 
ically chosen to compare to the results of iRahal et al.l 
IpOOa i. In this case, the constraints on the combination 
c*nm-c*bd are ~ ^ times weaker primarily because in this 
direction only disk lenses are contributing to the event 
rate. 

How well can we determine the minimum mass of a 
nomad, m m ; n , from a WFIRST type survey? Because a 
given mass nomad will produce events over a fixed range 
of timescales (for an assumed velocity distribution func- 
tion) the answer to this question depends on the value 
of TO m i n itself. If the mean event timescale at a given 
TO m in is significantly less than the cadence of the survey, 
then observations will not effectively be able to probe 
this mass scale. 

In Figure[5]we show the resulting fractional uncertainty 
on m m i n for a cadence of 15 minutes and bulge observa- 
tions. In this case, for m m - m > 10 _5 M Q we find frac- 
tional uncertainty > 30%. In fact down below the Earth 
mass scale for m m ; n > 10~ 6 M Q we still find fractional 
uncertainty ~ 50%, below which there is degradation of 
the constraints because the event rate in the observable 
timescale window becomes too low. 

6.2. Large-Scale Surveys: Gaia and LSST 
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Figure 4. Left Panel: Joint constraints on the slope of the mass function in the nomad region, ct nrn , and the slope in the low mass 
stella r regime, ctbd for b ulge sources. Unfilled contours assume the exposure and efficiency of modern bulge surveys, and match the results 
from[Sunii ct al. (2011]). Inner contour is 68% c.l. and outer contour is 95% c.l. Filled contours are the projected constraints for an 
exposure and efficiency expected for WFIRST. Right panel: Similar to filled contours on the left, except assuming disk sources. 
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Figure 5. Fractional uncertainty on m m i n for a cadence of f5 
minutes and bulge observations. 

We now extend to consider projected constraints on 
«nm-«bd from all-sky observations. In this case, estimat- 
ing the detection efficiency of the survey will be crucial 
in order to understand what fraction of the total event 
rate shown in Figure [3] will be accessible. 

The two primary templates we consider for large scale 
surveys are those being planned for Gaia and LSST. 
These multi-purpose surveys are not expected to have ca- 
dence as high as the dedicated inner Galaxy observations 
discussed above, so they will not be as sensitive to very 
short timescale microlensing events. However by their 
nature all-sky observations do probe the nomad popula- 
tion on a Galaxy-scale that are inaccessible to dedicated 
pointings towards a fixed region of the Galaxy. 

6.2.1. Gaia 



As our first example of a survey with a non-uniform 
cadence, we consider Gaia, which is scheduled to launch 
in 2013. Though Gaia is primarily designed as an astro- 
metric mission, it will have a single measurement pho- 
tometric accuracy of ~ 10 mmag for sources brighter 
than its broadband 20th magnitude. Because of the Gaia 
scanning strategy, the sampling for each star is not uni- 
form during the mission lifetime. Measurements will be 
grouped into epochs, during which an observation is per- 
formed in ~ 6 hr intervals. The mean number of mea- 
surements per epoch is ~ 5, though some epochs will 
have a minimum of 7 measurements ([Ever fe Mignardl 
[2051) . The mean number of visits between epochs is 
25-35 days, though depen ding on Galactic l atitud e we 
estimate from the results of lEver fe Mignardl ([20051 ) that 
~ 10% of the stars will have ~ 5 days between epochs. 

Motivated by these specifications, we model Gaia ob- 
servations by considering a quasi-irregular sampling pat- 
tern. For the baseline model we assume 25 days between 
epochs, and within each epoch there are five photometric 
measurement s. This is the approxim ate mean sampling 
rate of Gaia ([Ever fe Mignard|[2005l) . The survey is run 
for a total of T f, s = 5 years, resulting in a mean of 300- 
400 photometric samples for the lifetime of the survey. 
To model the distribution of disk sources we use the V- 
band luminosity function d escribed above, along with the 
IBelokurov fc Evans! ([2002D dust extinction parameteriza- 
tion. 

For a Gaia-likc sampling, the short timescale events, 
tE ~ 1 day, will occur during an epoch, and it is pos- 
sible that a peak of the microlensing event will not be 
discernible. To account for this, we modify the detection 
criteria. For a simulated event at an input timescale, we 
again search for three points on the lightcurve that are 
greater than 3-tr deviations from the baseline magnitude 
of the source. In addition we include a second, stricter 
cut to the detection criteria, namely that the peak of the 
event is observable. 

Given the above algorithm, for the Gaia sampling pat- 
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tern, at = 1 day we find a detection efficiency of 1%. 

For the Gaia cadence and estimated efficiency, in Fig- 
ure [6] we show the joint constraints on a nm and a^a for 
all-sky observations. Here we have assumed a five year 
lifetime of the mission. In this case the constraints are 
similar to the current constraints on these parameters 
because of the similar event rates after our detection ef- 
ficiency cuts have been implemented. 

6.2.2. LSST 

As our second example, we examine the somewhat 
deeper survey we anticipate being carried out by the 
LSST. This system will repeatedly survey the entire 
visible southern sky to a 5-sigma point source depth 
of g = 25.0, r = 24.7 per epoch. LSST is ex- 
pected to have a mean cadence (across all filters) of 
less than 4 days and a mission lifetim e of 10 years 
()LSST Science Collaborations et ai]|2009t) . To achieve a 
synoptic survey, the LSST will follow a logarithmic sam- 
pling pattern, with 15-second exposures separated by 15 
seconds, 30 minutes, 3-4 days and one year, with con- 
siderable scatter in the two intermediate cadences to al- 
low flexible scheduling. Two back-to-back exposures con- 
stitute a "visit"; the baseline plan has each field being 
visited twice on any of its observation nights. As with 
Gaia, detection of a nomad by microlensing requires see- 
ing both sides of a peak in the lightcurve, suggesting that 
events with timescales less than 3 days may be difficult 
to detect. 

To approximate the sampling strategy of LSST, we as- 
sume a uniform cadence for the lifetime of the survey. 
From the formalism above we have calculated the de- 
tection efficiency for cadences of both 1 and 4 days; we 
consider higher cadence dedicated campaigns with LSST 
below. We find that only for a 1 day cadence is it possible 
to achieve 1% detection efficiency for timescales is > 1 
day. For a 4-day cadence, the efficiency for detecting no- 
mads drastically drops (though in this case a large num- 
ber of brown dwarf events will still be measured very 
precisely). We utilize this optimistic 1% efficiency when 
we calculate the projected constraints on a nm -abd for the 
uniform cadence model. 

To model the distribution of sources for our LSST 
predictions we use the /-band luminosity function from 
IZheng et al.l (|200l . In this case, dust extinction mod- 
eledby adopted in the model of iBelokurov fc Evani 
(|2002D and scaling accor ding to the standard ext inction 
law between wavebands (|Rieke fc Lebofskvlll985l) . 

The left panel of Figure |6] shows the results of the anal- 
ysis. As indicated, the constraints are weaker relative to 
the left panel; this is mainly due to the reduced efficiency 
as compared to the Gaia sampling model. 

6.3. Kepler 

As our final example we study the nomad event rate 
measurable by the Kepler satellite Kepler monitors 
~ 100 deg 2 towards the Cygnus region, and has a pho- 
tometric precision of approximately 80 ppm for sources 
brighter than V = 13, and a few percent for sources at 
V = 20. Kepler is complete down to V = 17. The in- 
tegration time is 30 minutes for the majority of Kepler 



sources. For a n 



2 and m„ 



10 _8 M Q , and as- 



suming m m in = 10 down to 20th magnitude, the raw 
event rate is a few per year for 30 minutes < < 1 day. 
For a nm < 2 we find less than one event per year. How- 
ever, these predictions are for uj< = 1; due to its precise 
photometry of ~ 80 ppm for bright stars with V < 13 
the event rate is proportionally larger for ut > 1. Thus 
discovery of anomalous microlensing events in the Kepler 
data may indicate a steep value for the nomad mass func- 
tion, and warrant a dedicated analysis of the photometry 
of Kepler stars. 

7. DETECTING SHORT TIMESCALE EVENTS 

In the above analysis we have restricted only to events 
with timescales sufficiently long to be detectable accord- 
ing to the criteria above. What if we relax this criteria, 
and expand to consider events with shorter timescales, 
over which the lightcurves are much more sparsely sam- 
pled? Is it still possible to detect microlensing events 
from lighter nomads over much shorter timescales? 

As an example let us consider the planned survey strat- 
egy of LSST. In a given LSST filter, each visit will consist 
of two consecutive 15 second exposures separated by a 
4 second readout interval. When possible, each field will 
be observed twice, with visits separated by ~ 15—60 min- 
utes. For stars with r < 20, the single visit photometric 
precision of each measurement is ~ 10 mmag. Though 
per visit the photometry is very precise and the 30 sec- 
ond cadence is short, two points on a lightcurve are not 
adequate to claim the detection of a microlensing event. 
However, the lack of variation of a source star over this 
timcscale in between visits could allow us to bound the 
existence of nomads with characteristic timescales ~ 30 
seconds. Assuming both the lenses and the sources to be 
in the disk, this timescale corresponds to a lens of mass 
< 10- 10 M o . 

Will lenses with such a small mass cause noticeable 
brightness fluctuations in a star when accounting for 
finite source effects? For our uniform source surface 
brightness model, we find that this depends on the lens- 
ing geometry. For example, with Dl = 1 kpc and D$ = 5 



kpc, a lens of mass 
brightness of ^4f s 



10- y ,10- 10 ]M o will have a peak 
1.01,1.001]. For lenses nearer to 
the source, Af s is reduced from these values. Though 
these arc smaller brightness fluctuations than typically 
searched for in microlensing events, the presence of these 
objects may be limited given the photometric precision 
of LSST. 

In perhaps more near of a term, we may entertain the 
prospect of a dedicated telescope that is capable of de- 
tecting short timescale microlensing events that last for 
as little as tens of seconds. As an example, we will con- 
sider a liquid mirror telescope similar in design to the 
six-meter Large Zenith Telescope Q, only in our case po- 
sitioned in the Southern hemisphere to cover the Galac- 
tic bulge. For a 4000 x 4000 CCD chip with l" per pixel 
and a 6-m aperture, in less than a day a patch of area 
~ 100 square degrees could be scanned. We focus on the 
/-band, and take the bulge as an exa mple with an sur - 
face brightness of 17.6 mag arcsec" " 2 (|Terndrupl[l988[ ). 
Assuming that the signal-to-noise is dominated by the 
unresolved light from the bulge and shot noise, for a star 



6 http://kcplcr.nasa.gov/ 



http: / /www. astro. ubc.ca/lmt/lzt/ 
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Figure 6. Joint constraints on a nm and «t,d f° r all-sky observations. The left panel shows results for the irregular sampling pattern that 
represents the Gaia mission, while the right panel represents the uniform sampling pattern that may be achievable with LSST. 



with / = 19 the signal-to-noise is S/N ~ 37-\/i/sec. A 
10-second exposure then gives a photometric precision of 
~ 1%, and during this time a star crosses through ~ 10 3 
pixels. This would likely be sufficient to obtain several 
points on a lightcurve to measure a microlensing event 
with tE ~ 30. sec. 

For a nomad mass function of dN/dM ~ M~ 2 with 
m-min = 10 _9 M Q , for a 100 deg 2 patch that passes 
through the Galactic center we find ~ 50 events per year 
with ts > 30 seconds for source stars with I < 19. The 
event rate may even be up to an order of magnitude 
larger for steeper values of the mass function over the 
range 10 _9 M Q . It is also worthwhile to point out that 
a telescope designed along these lines could also be a 
relatively inexpensive endeavor. Further, a liquid mirror 
telescope with mercury could extend to the near-infrared, 
where the reflectivity would be similar to that in the op- 
tical. 

8. DISCUSSION AND CONCLUSIONS 

We have estimated that there may be up to about 10 5 
compact objects per main sequence star in the Galaxy 
that are greater than the mass of Pluto. A dedicated 
high cadence survey of the inner Galaxy, such as would 
be possible with WFIRST, could measure the number 
of nomads greater than the mass of Jupiter per main 
sequence star to ~ 13% , and the corresponding number 
greater than the mass of Mars to ~ 25%. Also WFIRST 
can measure the minimum mass of the nomad population 
to about 30%. Large-scale surveys, in particular that of 
Gaia, could identify nomads in the Galactic disk that are 
greater than about the mass of Jupiter. 

Observations along the lines that we discuss will con- 
strain the nomad population of the disk relative to 
the bulge, and will also more generally improve the 
star-star microlensing event rate in the disk and the 
solar n eighborhood, about which very l i ttle is now 
known ijGaudi et alJI^llFukui et al.ll2007t iRahal et al.l 
I2009D . Further, they will improve our understanding of 
the mass function of low mass brown dwarfs and supcr- 
Jupiters, and the distinction between these classes of ob- 



jects ([Spiegel et al.ll201lD . 

How will these measurements compare to modern mi- 
crolensing measurements of low mass brown dwarf pop- 
ulation from disk observations? To answ er this ques- 
tion w e can briefly consider the results from IRahal et all 
(|2009f ) . These authors find a total of ~ 20 events in three 
fields in which the lenses are primarily disk sources, and 
in particular there are two very short timescale lenses, 
at tE = 7 and 12 days. While the data are not suffi- 
cient at present to perform a full statistical analysis and 
constrain a nm and abd, from an analysis of these data 
one may deduce that a steeper model brown dwarf mass 
function is favored over a more shallow model. The in- 
clusion of the nomad population does mildly improve the 
statistical fit, though in order to probe this population 
with disk observations a survey must build up a sufficient 
event rate in the ~ 1 — 10 day timescale bin. 

How will the future microlensing measurements we 
discuss compare to dir ect measureme n ts of the brown 
dwarf mass function? iMetchev et~ai1 (|2007j ) find that 
for warm brown dwarfs the mass function may be flat, 
Oibd = 0. For cooler brown dwarfs the recent WISE ob- 
servations are consistent with a w ide range of «bd be- 
tween — 1 (jKirkpatrick et aL1l201lD . Other microlensing 
observations shed light on the brown dwarf mass func- 
tion , though they do no t clarify the picture. For exam- 
ple, iGould et all (f2009h uncover an extreme magnifica- 
tion microlensing event and interpr et it as due to a thic k 
disk brown dwarf. According to IGould et al.l (|2009fl . 
there is a very low probability to observe this event given 
our standard population of brown dwarfs in the Galactic 
disk given the large velocity of the event. The existence 
of these events either implies that we have been lucky 
to observe events at all (in particular with the large ob- 
served magnifications), or that the local population of 
low mass and low luminosity stellar remnants is larger 
than is presently predicted. 

If a nomad is identified via the methods described in 
this paper, there are a number of follow-up observations 
that are possible. For example even though Gaia will 
only do on average 1-2 one-dimensional astrometric mea- 
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surements per epoch, it may be possible to confirm the 
photometric detection with astrometry for the brightest 
sources by comparing the centroid of the source during 
the event to the baseline centroid as determined over sev- 
eral epochs during the course of the mission. 

Finally we note that an additional outcome of the ob- 
servational approach discussed above, especially regard- 
ing the detection of short timescale microlensing events, 
is that upper limits may be set on the density of nomads. 
This could set very interesting constraints on the popu- 
lation of planetesimals in nascent planetary systems. 
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